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Cushion  beam,  feet 
Appendage  drag  coefficient 
Skin  friction  coefficient 
Total  vehicle  drag,  lbs 
Appendage  drag,  lbs 
External  aerodynamic  drag,  lbs 
Dry  stores  duration,  days 
Seal  (or  ski)  drag,  lbs 
Sidewall  frictional  drag,  lbs 

Additional  (external)  sidewall  frictional  drag,  lbs 
Sidewall  wavemaking  drag,  lbs 
Wavemaking  drag,  lbs 
Froude  number  based  on  £ 

Wavemakir.g  drag  parameter 
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Acceleration  of  gravity,  ft/sec 
Fully  wetted  sidewall  depth,  feet 

Waterjet  system  loss  coefficient  (includes  internal  and 
elevation  losses  plus  an  equivalent  loss  to  account  for 
the  additive  inlet  induced  external  drag) 

Lift-drag  . atio 

Cushion  length,  feet 

Number  of  accommodations 

Design  cushion  pressure,  psf 

lhaft  horsepower  delivered  to  the  pump,  h.p. 

Power  required  at  . h.p. 
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Power  required  per  engine  at  VrK  , h.p. 
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Lift  horsepower  required,  h.p. 
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Available  maximum  continuous  shaft  horsepower,  h.p. 

Available  maximum  intermittent  shaft  horsepower,  h.p. 

Propulsive  horsepower  required,  h.p. 

1.13  times  power  required  at  , h.p. 

max 

Design  foil  loading,  psf 

Air  dynamic  pressure,  psf 

Water  dynamic  pressure,  psf 

Vehicle  range,  miles 

Wavemaking  resistance  coefficient 

2 
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Maximum  vehicle  speed,  smooth  water,  knots 

Take-off  speed,  knots 

Vehicle  gross  weight,  lbs 

Auxiliary  systems  and  outfit  and  furnishings  weight,  tons 
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Fuel  weight,  tons 

Propulsion  system  weight,  tons 

Payload  weight,  tons 

Hull  structure  weight,  tons 

Lift  systems  weight,  tons 
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Ratio  of  V to  V 
K * K 

d max 
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Lift  horsepower-to-propulsive  horsepower  ratio 
Overall  propulsive  coefficient 
Pump  efficiency 

Maximum  pump  efficiency 
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INTRODUCTION 

The  project  objective  is  to  reccunend  "Reference  Craft"  for 
use  in  an  exploratory  development  program  of  waterjet  propulsion 
systems.  The  vehicle  concepts  recommended  in  this  report  are 
limited  to  200  tons.  However,  other  vehicle  sizes  may  readily  be 
examined  by  utilizing  the  method  presented  herein. 

A rapid  method  is  developed  for  estimating  the  overall  propulsive 
coefficient  and  the  drag  of  a vehicle  concept  over  its  operating  speed 
range.  Thus,  an  estimate  of  the  shaft  horsepower  required  as  a function 
of  speed  may  be  obtained.  Preliminary  weight  estimating  relationships 
are  also  included  so  that  predictions  of  fuel  weight  and  payload  weight 
may  be  made.  A computer  program  utilizing  this  method  and  applicable  to 
SES  vehicle  concepts  is  contained  in  Appendix  A. 

Some  of  the  equations  used  in  the  method  are  presented  without 
derivation.  Future  reports  will  cover  the  development  of  these 
relationships . 


THE  OVERALL  PROPULSIVE  COEFFICIENT 
FOR  WATERJET  PROPELLED  VEHICLES 

The  overall  propulsive  coefficient  n is  defined  by, 

D(1.6878Vr) 
n “ 550  ' P 

where  D is  the  total  vehicle  drag,  exclusive  of  that  drag  induced  by  the 
inlet,  in  pounds,  VT.  is  the  speed  of  the  vehicle  in  knots,  and  P is  the 
shaft  horsepower  delivered  to  the  pump  at  the  vehicle  speed  V^.  For 
purposes  herein  we  assume  e transmission  efficiency  of  unity. 
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(or  parametric  analysis  of  waterjet  propelled  high  performance 
vehicle  concepts  operating  in  a patrol  type  mission,  i.e.,  cruising  at 
some  specified  speed  (e.y.,  20  knota)  and  having  a specified  maximum 
(dsah)  speed  (e.g.,  75  knots),  it  is  necessary  to  be  able  to  estimate 
the  powering  requirements  at  the  specified  cruising  speed  so  that  fuel 
weight  requirements  corresponding  to  the  specified  cruising  range  tnay  be 
determined.  When  the  desired  cruising  speed  coincides  with  hull-borne 
or  pre-hump  operation,  it  will  be  assumed  that  the  vehicle  concept  haB  a 
hull-borne  waterjet  propulsion  system,  in  addition  to  the  primary  (high 
spaed)  system,  which  can  he  characterized  by  supposing  that  a small  hull- 
borne  inlet /pump/engine  y:'em  is  Included  aa  part  of  the  propulsion 
system.  The  equation  (to  • developed)  fc  r n will,  in  addition  to 
providing  an  estimate  of  n for  the  primary  waterjet  system,  enable 
characterization  of  this  "secondary"  propulsion  system  efficiency  over 
the  hull-borne  speed  range. 


An  expression  for  the  overall  propulsive  coefficient,  as  giv^n 
in  reference  (1)  lEquation  6],  is 


n/np  - 2 ( (l  - k - ^ )1/2  - i]v2/(2gH)  [i] 

v V 


where 

rip  is  the  pump  efficiency 

k is  the  ratio  of  total  inlet  and  internal  energy  losses 

(exclusive  of  elevation  changes)  to  the  vehicle  velocity 
head  (V2/ 2g) 

h is  the  difference  in  elevation  between  the  nozzle  exit  and 
the  water  surface,  ft 

H is  the  head  produced  by  the  pump,  ft-lb/lb 

2 

g is  the  acceleration  due  to  gravity,  ft/sec 

V is  the  forward  speed  of  the  vehicle,  ft/sec 
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As  no  attempt  la  made  herein  to  estimate  the  Inlet  Induced  external 
dreg,  we  account  for  this  gdditlva  drag  by  penalizing  the  ovetall 
propulsive  coefficient.  This  Is  done  by  replacing  K - k + 2*h  (the 
magnitude  of  which,  In  reference  (■)),  Is  taken  to  be  fixed  fi>r  a fixed 
geometry)  by  e waterjet  system  loss  coefficient  , which  Includes 
internal  (Inlet  plus  duct)  and  elevation  losses  plus  an  equivalent  loss 
to  account  for  the  additive  Inlet  induced  external  drag.  Equation  [i] 
then  becomes 


n/n. 


((1  ••  K.  ■«  )1/2  - l]V2/(gH) 


[2] 


As  pointed  out  in  reference  (1),  this  equation  has  a maximum  value  for 
a given  value  of  Kd  so  that  there  is  an  optimum  value  for  H.  iVs  is 
found,  as  in  reference  , by  differentiating  Equation  [2]  with  respect 
to  H*  “ 2am.  and  solving  for  H.  The  result  is 

\t ^ 
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opt 


T |Ki  + 


where  V denotes  the  particular  value  of  V for  which  H is  optimum, 
d 

Reference  (1)  suggests  that  if  the  magnitudes  of  pump  head  and  flow  rate 
at  cruise  (design)  speed  are  maintained  at  take-off  speed,  then  the 
additional  thrust  required  at  take-off  may  be  produced.  Following  this 

reasoning,  we  assume  that  the  magnitude  of  H is  a constant  H - Hopt’ 

2 P 
V, 

H - — (K.  + ylT]  if  0 < V < V 

g d a — a 


Replacing  H in  Equation  (2]  by  this  relationship  and  setting  V « aV 
r a max 

0 < a < 1,  gives 


n/nT 


(«  - Kd  ♦ 


2<*d  + 

<Vav 

max 


1/2  , jV“VK  > 

J - 1 ] max 

<Kd  + 


[3] 


if  0 < V 1 V 


where  the  vehicle  speeds  are  now  given  in  units  of  knots. 


An  estimate  for  the  pump  efficiency  is  needed  which  has  a 
specified  maximum  np  - np  when  . Moreover,  the  effect  of 

should  be  included  since  fffe  magnitude  of^the  total  head  above  vapor 
pressure  at  the  impeller,  Hgv  (referred  to  ae  the  net  positive  auction 
head),  is  given  by 


SV 


(1  - k) 


2g 


+ h 


SV 


Pi 


where  is  the  difference  between  atmospheric  and  vapor  pressure  expressed 
as  a head  in  feet  of  water  and  h is  the  elevation  of  the  pump  inlet 
acove  the  free  water  surface  (' f the  pump  inlet  is  below  the  free  water 
surface,  then  hp,  is  a negative  quantity).  Referring  to  the  equation  for 
Hsv  it  is  seen  that,  for  a fixed  vehicle  speed,  the  magnitude  of 
increases  as  k (or  K^)  decreases  with  a corresponding  increase  in  the 
pump  specific  speed  for  a given  pump. 


Test  results  for  the  XR1-B  and  calculations  were  used  to  determine 
the  variation  of  pump  efficiency  with  vehicle  speed.  This  variation 
can  be  characterised  by  a slowly  increasing  function  of  vehicle  speed 
which  reaches  a maximum  value  at  design  speed.  Such  a function  is 

Vk  2 

np  - np  (l  - c (l  - n 

max  K. 

a 

where  np  is  the  maximum  pump  efficiency  and  is  assumed  to  be  a 

function”??  K. . 
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We  assume  that  the  above  functional  relationship  between  pump 
efficiency  and  vehicle  speed  is  generally  true  and  that,  as  discussed 
previously,  the  propulsion  system  losses  (as  reflected  by  K^)  slightly 
affect  the  pump  efficiency.  To  estimate  this  effect,  the  ratio 
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r*  rw jir  r*  tb« 


v^/(l  +•  »^fj}  la  asaumed  for  the  constant  as  it  provides  ■'.nail 

variations  in  the  magnitude  of  nD  for  a considerable  range  of  K . values. 

* a 

Th«  validity  of  thio  assumption  will  have  to  be  verified  when  more 
experience  is  gained  in  the  application  of  water jet  propulsion  systems 
to  high-speed  vehicles.  The  equation  selected  to  estimate  the  pump 
efficiency  is 


[i  ‘ 1 + (1  ' aVv. 


y ) a o < vK  < vK 


Combining  Equations  [3]  and  [4]  gives  an  estimate  of  the  nveiui.1 
propulsive  coefficient  p over  the  speed  regime  (0  < V < V ) for  a 
vehicle  concept  with  a specified  loss  level  (K^)  and  a spec¥ffed  maximum 


pump  efficiency,  n 


. i.e. , 


if  0 < vK  < vK 


For  high-speed  operation,  the  design  speed  -\r  r " aV  , should  be  set 

‘Snax 

equal  to  V , i.e.,  a ■ 1.  If  a secondary  propulsion  system  is 
Nnax 

assvT .d  for  purposes  of  cruise  range  specification  at  low  speed,  then 
a second  set  of  values,  a.  K^.  and  n,,  . . should  be  selected  and 

used  in  Equation  (5]  to  determine  themnXfor  this  propulsion  system. 


In  using  this  method  for  estimating  the  overall  propulsive 

coefficient  of  waterjet  propelled  vehicle  concepts,  the  water jet 

system  loss  coefficient  K may  be  treated  as  a parameter. 

a 


In  this  case  the  following  tentative  range  of  magnitudes  for  K 

d 

are  recommended  for  preliminary  performance  estimates: 


Water  let  System  Reconmended  Values 

Flush  or  semi-flush  inlet  (hull)  .4  - .5 

Hydrofoil  strut/pod  ram  inlet  .65  - .75 

Semi-f lush/ram  inlet  (short  strut)  .5  - .65 

SES  DRAG  AND  ®OWER  ESTIMATES 

The  drag  estimates  for  SES  (rigid  sidewall  ACV's)  are  obtained  by 

a summation  of  component  drag  relationships  nased  essentially  on  the 

procedure  given  in  references  (2)  and  (3),  Aa  an  approximation,  it  is 

assumed  that  the  SES  is  fully  cushion  supported  (i.e.,  W » p S,  where  W 

o 

is  the  gross  weight  in  lbs,  p is  the  design  cushion  pressure  in  psf, 

0 2 

and  S i?  the  cushion  clan form  area  in  ft  that  the  cushion  piantorm 
1b  rectangular,  that  the  sidewall  length  is  equal  to  the  cushion  length 
(H) , and  that  the  daylight  clearance  is  zero.  Consistent  with  the  latter 
assumption,  the  ram  (or  momentum)  d.:^  is  taken  to  be  zero. 

The  equations  used  to  predict  the  component  drag  to  weight  ratios 
are  09  follows: 

Wavemaking  drag,  D^:  From  reference  (2)  , this  component  is  given  by 


.)  , p 

_W  . ± . _£L  . f [6] 

Pg  l i 


where  p / g,  is  the  design  pressure  to  cushion  length  ratio  and  f is  the 
v x 

wave  drag  parameter.  A cloaed  form  analytical  expression  for  f„  has 

X 

been  developed,  assuming  a uniform  pressure  distribution  acting  over 
a rectangular  region.  It  Is 
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where 

A - fc/b  Is  the  le-igth-beam  ratio  cl  the  rectangular  planlortn 

y - .577...  is  Euler's  constant 
F^  « V/r^gl  is  the  Froude  number 

F,H  ■ V^ST  in  the  primary  hump  Froude  number , which  is  estimated 


by  F 


01745  A2 


Sidewall  frictional  dreg,  D 

#,« 


For  2 sidewalls,  this  component  is  given  in  reference  (?)  by 


where  C,  is  the  kin  friction  coefficient,  q is  the  water  dynamic 
t w 

pressure  (2.P35V^),  and  h is  the  fully  wetted  sidewall  depth.  The 

IN  ft 

following  estimate  for  h is  used. 
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s,v 

reference  (2),  this  component  la  given  as 
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Sidewall  wavemaking  drag, 


D : 
s,w 


Thl9  component,  as  Riven  in  reference  (1)  , la 


r>  2 2 

L9,W  B f 9) 

W " W ~ir  L 

Values  of  the  wavemaking  resistance  coefficient,  r,  for  a full  form,  having 
B/L  - .C265,  B/H  - .183,  and  L/H  * 7,  were  taken  from  reference  (4).  For 
Froude  numbers  between  0 and  1.80,  the  magnitude  of  r is  computed  by  inter- 
polation of  these  stored  values.  For  larger  Froude  numbers,  the  following 
approximation  is  used. 

r ■=  6.125  • ~^2  ‘ ln  F£  lf  1'80  I F£ 

Ff 

The  assumption  is  also  made  that  H ■ h . 

a 

External  aerodynamic  drag,  D^:  As  given  in  reference  (2)  , this  component 

is 


D 

e 

w 


P 
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[10] 


where  q^  is  the  air  dynamic  pressure  < .00339VK> . Here  we  have  assumed 
Che  estimate  of  .2/(i/b)  for  the  aerodynamic  drag  coefficient  as  given  in 
reference  (5)  . 


Sesl  (or  ski)  drag,  D : For  this  component,  a flat  planing  surface  is 

assumed.  The  average  bottom  planing  speed  is  taken  to  be  equal  to  the 
vehicle  speed  and  the  ski  width  is  assumed  equal  to  b.  The  drag  for  two 
seals  is  then 
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•1 


DSK-2  • 2*V  Cf£,'  £'b  " 2^  C^t’b 


where  Cf  is  the  flat  plate  akin  friction  coefficient  and  t'  la  the 

chord  length  of  the  eki.  Assuming  *'  • .0251  and  Cf  ■ 2Cf  gives 

‘l'  1 


. , . c . s 

W P 
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Appendage  drag,  D^:  Reference  gives  an  estimate  for  the  drag 

coefficient  of  parabolic  cross  section  appendages.  If  the  values  .01 
and  .1  are  assumed  for  appendage  area  to  cushion  area  ratio  and  thickness 
to  chord  ratio,  respectively,  then  the  drag  coefficient  given  therein 
becomes 


CnA  - .0000393  + .02C, 
DA  f 


The  appendage  drag  is  then  taken  to  be 


Da  CPAqw 


The  tocal  drag  to  weight  ratio  (D/W)  is  obtained  as  the  sum  of 
Equations  [6]  through  [12J,  i.e.. 


DfI  D D D D h „ D 

D m _W  s,a  s ,v  s,w  e SK  a 

w " w w w * w w w w 


from  which  we  obtain  the  propulsive  horsepower  to  weight  ratio  as, 


P _ V„ 

p m D K 

W " W 326n 


where  is  the  vehicle  speed  in  knots  and  n is  the  overall  propulsive 
coefficient . 


A 


It  la  assumed  chat  the  lift  horsepower  required  (P^)  le  proportional 
to  propuleivj  power  required,  i.e. 


(15J 


The  smooth  water  shaft  horsepower  required  is  now  given  as 


P - 


(16) 


where  W is  the  gross  weight  of  the  SES  in  pounds. 


HYDROFOIL  DRAG  AND  POWER  ESTIMATES 


\ 


Smooth  water  drag  estimates  for  hydrofoil  vehicles  are  based  on 
empirical  relationships  developed  for  estimating  vehicle  lift-drag  ratios, 
L/D .over  the  vehicle's  speed  regime  for  hydrofoils  having  either 
subcavitating  or  supercavitating  lift  systems. 

For  hydrofoil  vehicles  having  subcavitating  lift  systems,  the  L/D 

curve  (as  a function  of  vehicle  speed)  decreases  rapidly  with  increasing 

speed  in  the  nullbome  condition  until  take-off  (where  L/D  is  a minimum) 

after  which  the  L/D  curve  increases  to  a maximum  value  (depending  on 

cavitation  conditions)  and  then  decreases  again  as  cavitation  becomes 

more  severe.  To  include  the  effect  of  various  take-off  speeds  on  L/D, 

the  take-off  speed,  V is  defined  by  setting  VK  ■ BVv  ,0  <6  <1. 

*t  max 

To  provide  the  variation  of  L/D  over  the  whole  vehicle  speed  range, 
the  following  variation  for  subcavitating  L/D' a is  selected 
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Mil ■ ijjfc. ■■  JfcJ fa  liUMCA.. | r I ill1  mi  -I, tfn 


550 


1 


Z' 


,2.5 


max 


L/D  -< 


550  1 


\ e1-5 

max 


550 


<vvk  y 

max 


(-JL_)2 

'v  ’ 

K 

max 


V..  <V„/V„  ) 

max 


V K N'K' 'K 

V max 


lf  0 < vK  < vK 


if  VKt  - VK  < ^ VK  (17] 


max 


lf  S VK  <-\<-  VK 

max  max 


where  Che  empirical  constant  Is  selected  so  that  the  L/D  at  50  knots  Is 
11.  This  Is  based  on  an  average  of  L/D’s  attained  by  existing  hydrofoils. 

The  variation  of  L/D  with  vehicle  speed  for  hydrofoil  vehicles 
with  supercavitating  lift  systems  Is  assumed  to  be  similar  to  the  variation 
of  L/D  with  vehicle  speed  for  subcavitatlng  lift  systems  throughout  the 
bullbome  speed  range.  After  reaching  a minimum  at  take-off,  the  L/D  is 
assumed  to  slowly  Increase  throughout  the  foilbome  speed  range.  The 
following  variation  of  L/D  with  vehicle  speed  reflects  the  above  assumptions. 


. s2-5 


L/D 


max 

530 

max 


<vK/vK  )- 

max 


<VVK  > 

max 


1/2 


if  o < VR  < V 


if  V < V < V 
K*.  - K - K 
t max 


[18] 


The  same  empirical  constant  is  selected  as  It  gives  L/D’s  of  7.33  and 
5.5  at  75  and  100  knots,  respectively.  These  magnitudes  are  considered 
to  be  attainable  with  careful  design  in  the  near  future. 


The  smooth  water  shaft  horsepower  delivered  to  the  pump  is  given 

as , 


6.87  t VK 
n (L/D) 


(19) 


11 


where  A - W/2240  ia  the  gross  weight  of  the  hydrofoil  in  tons. 

PLANING  CRAFT  DRAG  AND  POWER  ESTIMATES 

The  smooth  water  drag  estimates  for  planing  craft  are  based  on 
Series  62  data  (reference  (7^1.  For  the  planing  craft  considered 
herein.  Figure  12(a)  of  reference  (7)  was  used  to  determine  the 
resistance  to  weight  ratios  D/W,  i.e.,  a length  to  beam  ratio  of  5.5 
was  selected  and  the  82  LCC  curves  were  used. 


The  shaft  horsepower  required  is  given  by 


P 


D 

W 


n 


120] 


WEIGHT  ESTIMATES,  ALL  VEHICLES 

In  order  to  obtain  first  order  payload  weight  estimates,  it  is 
necessary  to  estimate  the  component  weights  o'"  the  selected  vehicle 
concepts.  Ordinarily,  one  would  specify  payload  weight  and  estimate 
the  grobs  weight;  however,  in  this  case  an  iterative  procedure  19 
required  to  solve  the  weight  equation  for  gross  weight.  By  specifying 
a gross  weight  (A)  the  payload  weight  may  be  calculated  explicitly. 

The  following  relationships  were  selected  for  use  In  estimating 
component  weights. 

Hull  Structure;  The  hull  structural  weight  group  W includes  the 

" 1 8 

main  body  and  hull  structure,  superstructure,  and  machinery  foundations. 
A minimum  structural  weight  envelope  given  in  reference  was  selected 
for  hydrofoils  based  on  the  best  currently  available  material  technology 
and  construction  techniques.  This  was  increased  by  approximately  282 
for  planing  craft  based  on  a recent  planing  craft  design,  A weight 

estimate  for  SES,  as  a function  of  l/b  and  P/l,  was  developed  based  on 

o 

the  variations  presented  In  both  references  (6)  and  (9)  . 
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The  hull  structural  weight*  are  given  by 


.15  4+  .25  4’ 


W$  . 19  4 + . 32  4* 


; HYDROFOILS 


; PLANING  CRAFT 


"fr  l~F~> + 77 7175 


} ; SES 


; TONS 


Lift  Systems:  The  lift  systems  weight  group  Wa  , includes  the  foil 

system  weight  in  the  case  of  hydrofoils  and  the  seal  system  weight  in 
the  case  of  SES's.  The  equation  selected  for  estimating  the  weight  of 
hydrofoil  lift  systems,  whose  form  is  taken  from  reference  (6).  is 

assumed  to  depend  on  gross  weight  and  the  design  foil  loading  P , In 

2 8t 
lbs/ft  . The  lift  system  weight  for  SES's  is  taken  proportional  to 

gross  weight. 


W 


1^-  + jpll  a3/2  ; HYDROFOILS 
sf  *sf 

.054;  SES 

0;  PLANING  CRAFT 


; TONS 


The  magnitude  selected  for  Pgf  should  not  exceed  2100  PSF  as  the  resulting 
weight  estimate  might  be  too  low. 


Propulsion  System:  This  weight  group  includes  the  prime  movers 

and  the  transmission  and  propulsor  systems.  Basec  on  available  data 
for  two  recent  high-performance  vehicles,  SES  100A  and  PHM,  a specific 
weight  (lbs  of  propulsion  system/maximum  intermittent  shaft  horsepower) 
of  3 Ib/h.p.  la  selected  to  provide  an  estimate  of  marine  gas 
turbine/water Jet  propulsor/light  weight  transmission  propulsion  systems. 
If  a secondary  propulsion  system  is  required,  then  an  additional 
.5  lb/h.p.  is  assumed.  This  gives 
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w 


l 


001339  P^;  it  1 propulsion  system 


•001S62  P^ ; if  2 propulsion  systems 


; TONS  [23] 


where  PMI  is  the  available  maximum  intermittent  shaft  horsepower  (80°F) , 

The  magnitude  of  PMI  is  selected  such  tnat  it  exceeds  the  larger  of 

either  1.15  times  the  power  required  at  V,.  , (PgH  - 1.15P(VK  )),  or 

max  max 

the  power  required  at  hump,  PfV^)  or  take-off,  P(VK  )•  If  existing 
(or  planned)  gas  turbines  are  not  considered,  then  the  estimate 


may  be  used. 


PMI  • P<W] 


Electric  System:  This  group  is  a small  percentage  of  gross  weight, 

and  for  our  purposes,  the  following  equation  is  selected  based  on  an 
approximate  average  of  the  SES  100A  and  PGH-2 . 


w 


.03  TONS 


[24] 


Auxiliary  Systems  and  Outfit  and  Furnishings:  This  group  comprises 

a larger  percentage  of  gross  weight  than  the  electtlc  plant,  however, 
it  is  still  comparatively  small  relative  to  the  previous  groups.  The 
weight  of  thi6  group  is  taken  to  be  10  percent  of  the  gross  weight  on 
the  basis  of  the  SES  100A  and  a recent  planing  craft  design,  i.e.. 


W - .1A;  TONS 
aux 


[25] 


Complement,  Personal  Effects,  and  Stores:  This  group  consists 

of  the  accommodations  dependent  weight  items,  some  of  which  may  be  a 
function  of  duration  at  sea  measured  in  terms  of  dry  stores  duration 
in  days.  The  complement  and  personal  effects  weight  is  estimated  by 
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assuming  225  lbs/oan.  The  stores  weight  per  man  is  assumed  to  be 
proportional  to  the  dry  stores  duration.  Based  on  the  PGH-2  and 
a recent  planing  craft  design,  wa  obtain  the  following  equation 


W - .IN  + .015N  D ; TONS, 

C Co 


[26] 


where  N is  the  number  of  accommodations  and  D is  the  dry  stores 
c s 

duration  in  days. 

Fuel:  This  weight  group  includes  the  fuel  required  to  meet  the 

range  requirements  specified  for  the  mission.  The  fuel  weight  equation 
used  is 


_(sfc)R?e 

Uf  - A(1  - e2240£iV) ; TONS,  [27] 

where  sfc  is  the  specific  fuel  consumption  in  lbs/h.p.-hr.  and  R is 
the  specified  range  in  miles.  Estimates  for  sfc  are  made  by  using  an 
equation  developed  by  H.  D.  Marron,  Naval  Ship  Research  and  Development 
Center,  NSRDC,  Annapolis,  for  s'c's  corresponding  to  maximum  continuous 
power  ratings  of  second  generation  gas  turbines;  however,  as  an  approxi- 
mation we  assume  it  valid  for  the  power  per  engine  P where  ® is  the 

en  en 

estimated  power  required  at  V , P , divided  by  the  number  of  engines  in- 

K e 

e -.15 

stalled.  The  equation  is  sfc  = 1.9  P * . For  hydrofoil  and  SES  con- 

en 

cepts  where  hullborne  range  is  specified,  an  sfc  of  .4  is  assumed  to 
reflect  a smell  diesel  installation. 

The  payload  Wp/L  is  defined  to  consist  of  those  items  comprising 
NAVSHIPS  weight  groups  4 and  7 (communications  and  control  and  armament, 
respectively)  all  variable  load  payload  and  margins.  The  payload  weight, 
given  in  terms  of  Equations  [21]  through  [27]  via  the  weight  equation,  is 
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F/L  a 81  m e aux  ces  V * 


APPLICATION  OF  THE  METHOD  OF  ESTIMATING  THE  PRIMARY 
CHARACTERISTICS  OF  SES’S,  HYDROFOILS,  AND  PLANING  CRAFT 

The  equations  presented  for  estimating  the  powering  characteristics 
and  weight  characteristics  of  SES,  hydrofoils,  and  planing  craft  may  now 
be  solved  provided  certain  design  parameters  are  knewn.  In  this  section 
these  parameters  are  discussed  and  defined. 


Design  Specifications  Common  to  All  Craft 


These  are  essentially  selected  according  to  the  mission 
requirements.  They  are: 

V,,  Maximum  speed,  knots 

Max  Design  ..peed  (speed  at  which  maximum  propulsive  coefficient 

■y  J £ 

is  desired),  knots 

V Endurance  speed  (speed  corresponding  to  range  specification), 

knots 

R Range,  trlleb 

N Number  of  accomnodetions 

C 

D Dry  stores  duration,  day* 

8 

Prylrad  weight,  tone.  (.In  general  this  should  be  a 

specific ition  resulting  from  mission  requirements,  in 
which  case  Equation  [28]  becomes  an  implicit  equation 
In  gros.i  weight,  A which  can  only  be  solved  by 
Iteration.  For  purposes  herein  we  chose  to  specify 
A and  compute  for  eas a :.n  computation.) 


*Two  values  are  required  for  this  specification  If  the  endurance  speed 
specification  is  less  than  hump  or  take-off  speed. 


If 


SKS  Design  Parameters 


following  parameters  are  required  for  determining  characteristics 


of  SES: 
Symbol 


Specification 
groaa  weight 

cushion  length  to  beam  ratio 


square  of  maximum  Froude  no. 


Units 


non-dim 


design  pressure  to  length  ratio  non-dim 


non -dim  F„ 


Procedure 
for  Determining 

specify 

specify 

P /«.  - 2//E7b 

lit  2 

. 00678 (PQ/i)  VR* 


cushion  length 


design  cushion  pressure 


feet  £•  .08854  (VK  /F.  ) 


lba/f  t' 


K > 

max  max 


P - (P  /J.)  • t 
o o 


ratio  of  VR  to  VR 

d max 


non-dim 


waterjet  system  loss  coefficient  non-dim 


specify 


specify 


maximum  pump  efficiency 


non-dim 


■specify 


With  the  above  now  known,  Equations  [6]  through  [13]  are  used  to 
determine  the  drag/weight  ratio  at  any  desired  Froude  number;  0 < F^ 

Next  corresponding  magnitudes  of  VR/VR  are  determined  via 

sax 

v ^ 

V„  1 . 6878V„ 


and  then  Equations  [14]  through  [16]  are  used  to  determine  the  smooch 

water  shaft  horsepower  required  at  each  speed,  0 < V„  < V„  . The 

K Nnax 

magnitudes  of  the  component  weight  groups  ai • then  found  using 
Equations  [21]  through  [27].  The  payload  weight  then  follows  from 
Equation  [28],  If  the  payload  weight  is  insufficient,  then  the 
procedure  should  be  repeated  using  a new  gross  weight. 


If  secondary  propulsion  system  is  assumed,  then  two  values  may  be  required. 
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Hydrofoil  Design  Pirmltri 


The  following  parameters  ere  required  for  determining  the 
characteristics  of  hydrofoils: 


Procedure 

Symbol 

Specification 

Units 

for  Determining 

£ 

gross  weight 

tons 

specify 

£ 

ratio  of  to  V 

Kt  Kmax 

non-dim 

•y 

specify  (V  < 35) 
Kt  ~ 

P.f 

design  foil  loading 

* 

lbB/ft 

recommend 

1200  if  0 < V„ 

inax 

P , - 1600  if  55  < V.. 

8f  Snax 

2100  if  70  < VR 

max 

a 

ratio  of  to 

d max 

non-dim 

specify 

Kd 

wuterjet  system  loss 
coefficient* 

non-din 

specify 

max 

★ 

maximum  pump  efficiency 

non-dim 

specify 

If  0 < V < 55,  then  subcavitatlng  lift  systems  should  be  assumed 
end  L/D  estlma^fs  for  the  desired  speeds,  0 < £ 55  may  be  computed 

using  Equation  [17],  For  these  cases,  it  Is  recommended  that  £ - .5,  i.e., 
assume  that  the  take-off  speed  is  1/2  of  maximum  speed. 

If  V*.  > 55,  then  Equation  [18]  should  be  used  to  reflect  the 

max 

lower  l/D's  obtained  with  supercavltating  lift  systems.  If  < 70, 

max 

then  take-off  speed  may  be  assumed  to  be  1/2  of  maximum  speed,  i.e., 

6 * .5;  however,  for  70,  one  should  select  V - 35  and  determine 

, _ ,max  _ r 

the  appropriate  value  for  £. 

If  secondary  propulsion  system  16  assumed,  then  two  values  may  be  required. 
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Smooth  water  shaft  horsepower  magnitudes  for  each  are  given 
by  Equation  [19],  and  the  payload  weight  is  determined  by  use  of 
Equations  (21)  through  (28]. 

Planing  Craft  Design  Parameters 


The  parameters  needed  for  determining  the  characteristics  of 
planing  craft  are: 


Symbol 

Specification 

Units 

Procedure 
for  Determining 

A 

gross  weight 

tons 

specify 

LP/BPX 

length-beam  ratio 

non-dim 

specify 

s 

projected  chine  length 

ft. 

specify 

BPX 

maximum  beam  at  chine 

ft. 

specify 

K?/os^)2n 

bottom  loading  coefficient 

non-dim 

specify 

ed 

deadrlse  angle 

degrees 

12.5 

a 

ratio  of  to  Vr 

a max 

non-dim 

specify 

Kd 

waterjet  system  loss  coef*- 
f icient 

non-dim 

specify 

npmax 

maximum  pump  efficiency 

non-dim 

specify 

Magnitudes  of  resistance  to  weight  ratios,  D/W  should  be  obtained 
from  reference  (7)  after  having  selected  an  LCG  location  given  therein. 

The  procedure  may  then  be  completed  by  use  of  Equations  (20]  through  (28]. 
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CHARACTERISTICS  OF  WATERJET  PROPELLED 
HIGH  PERFORMANCE  VEHICLES 


Vehicle  concepts  considered  for  purposes  of  providing  s reference 
vehicle  for  the  selection  of  specif icetlona  for  a developmental  water Jet 
propulsion  system  are  categorized  according  to  a desired  maximum  smooth 
water  speed  capability  as  follovs: 


Vehicle  Concept 
Planing  Craft 
Hydrofoil 

v 5E= 

Hydrofoil 


Design  parameters  and  primary  characteristics  for  the  planing 

craft  concepts  are  shown  in  Table  1.  The  selected  design  parameters 

have  previously  been  definad.  Two  particular  sets  of  design 

specifications,  both  of  which  require  50  knot  maximum  smooth  water 

speed  capability,  are  selected.  One  is  for  a vehicle  concept  having  a 

50  knot  design  speed  VK  and  a 500-mile  range  R at  that  speed  V ; 

d *Se 

while  the  other  reflects  a concept  designed  to  be  most  efficient  at 

25  knots  and  to  have  a 500-mile  range  at  25  knots.  For  this  patrol 

mission,  however,  only  one  propulsion  system  is  assumed.  In  each  case, 

the  number  of  accommodations  N and  the  dry  stores  duration  D are 

c a 

assumed  to  be  12  and  7,  respectively. 


For  each  set,  gross  weights  A of  50,  100,  150,  and  200  tons  are 

selected.  For  reference,  the  planing  craft  PTF  has  a gross  weight  of 

approximately  85  tons  and  the  displacement  ship  PG  has  a gross  weight 

of  about  225  ^ons.  The  geometry  selected  is  the  same  for  both  sets. 

The  loss  coefficient  specified  for  each  set  is  .5,  which  is  assumed 

to  be  representative  of  waterjet  systems  having  flush  (or  semi-flush) 

Inlets.  Maximum  pump  efficiencies  n_  of  90%  are  assumed.  The 

rmax 
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TABLE  1 - 50  KNOT  PLANING  CRAFT  CONCEPTS  - DESIGN  PARAMETERS  AND  CHARACTERISTICS  (CONT.) 
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resultant  overall  propulsive  coefficient  n at  maximum  speed  Is  .527, 
while  at  25  knots  It  Is  .365. 


The  estimated  characteristics,  powering  and  weights,  are  presented 
next.  The  shaft  horsepower  required  by  the  pump  is  given  at  two  speeds, 


i.e.,  P(VK  ) and  P(Vj 


The  magnitude  of  PgH  is  also  given. 


The  number  and  model  of  possible  gas  turbine  Installations  along 

with  the  total  maximum  intermittent  PWT  and  maximum  continuous 

Ml  MC 

power  rating  (at  80®F  unless  otherwise  noted)  are  shown  for  each  concept. 


Tables  2 and  3 present  selected  design  parameters  and  resultant 

characteristics  for  the  50  and  75  knot  hydrofoil  concepts,  respectively. 

In  each  case,  two  "mission  profiles"  are  assumed;  the  first  specifies 

that  V„  - V and  that  the  range  at  Vv  m Vv  be  500  miles,  while 
Kd  Ke  max 

the  second  represents  a hull-borne  cruise  capability  of  2750  miles  for 

the  50  knot  hydrofoils  (2000  miles  for  the  75  knot  hydrofoils)  at 

V - . 2VK  " VK  with  a dash  capability  of  - For  this  patrol 

^d  max  e max 

mission,  a secondary  propulsion  pystem  is  assumeaT*  Parameter  values 

characterizing  It  are  given  in  parentheses  next  to  the  selected  parameters 

for  the  primary  propulsion  system.  For  all  cases,  the  number  of 

accommodations  N and  the  dry  store  duration  D are  taken  to  be  12  and  7, 

C 8 

respectively . 

For  each  set,  gross  weights  A of  50,  100,  150,  and  200  tons  are 

selected.  The  PCH-2  hae  a gross  weight  of  58  tons  and  the  PHM  is 

approximately  224  tons  gross  weight.  The  design  loading  F { selected  for 

the  50  knot  subcavitating  hydrofoil  concepts  is  1200  psf,  while  for 

the  75  knot  supercavitating  hydrofoil  concepts  it  is  2100  psf.  The  loss 

coefficient  K.,  taken  to  be  .75  for  the  primary  water jet  system,  is 
a 


A small  diesel  installation  is  assumed  but  not  identified. 
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TABLE  2-50  KNOT  HYDROFOIL  CONCEPTS  - DESIGN  PARAMETERS  AND  CHARACTERISTICS 
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o 

NO 

o 

rA 

CM 

»A 

<A 

O 

o 

US 

f** 

00 

•o 

CO 

tA 

Os 

• 

• 

• 

• 

• 

• 

'O 

»A 

Os 

00 

o 

NO 

CO 

NO 

© 

CM 

*4 

NO 

cn 

N 

*4 

CA 

N 

CM 

CM 

-4 

CM 

o 


m 

o 

*4 


O' 


o 

nO 


© 

<n 


m 

cm 


ua 


in 

*4 


o 

m 


m 


»A 

ao 


vO 

m 


O 

00 

nA 

NO 

<n 

Os 

r- 

o 

CM 

nO 

• 

• 

OO 

CM 

*4 

■»4 

rH 

*4 

rH 

ca 

NO 

NO 

r-. 

CM 

H 

rH 

o 

© 


ia 

N 


CA 

N 


*■4 

CM 


o 

*A 

00 

*4 

Mr 

CA 

ao 

SA 

o 

«A 

CM 

CA 

*— < 

<a 

co 

<M 

• 

- 

. 

• 

• 

• 

*-h 

CM 

SO 

NO 

m 

CM 

NO 

*— 1 

CA 

CA 

(A 

rH 

rH 

> 

A* 


TJ 

* 


K 

a 


a 


o 

to 


o 

CM 

O' 


SO 

fA 


O 

Os  « 

Os  O 


I 

u- 


o 

CM 


o 

o 

8 


CM 

•o 


>4 

I 


o 

NO 


o 

m 


o 

o 

oo 

NO 


*4 

I 


o 

m 


o 

CA 

s 


o 

o 

*4 


o 

SA 

lA 

nO 

ca 

O' 

o 

o 

o 

»A 

rH 

<N 

2 

o 

o 

CM 

fs 

r*«* 

*4 

• 

• 

• 

• 

• 

• 

• 

• 

i 

NO 

o 

CO 

*4 

•4 

rH 

o 

CA 

o 

CM 

tA 

CA 

H 

-4 

ao 

nO 

NO 

nO 

CM 

H 

rH 

H 

rH 

CM 

o 

f- 

NO 

KO 

cm 


o 

CO 

CO 

-4 

<*4 

<A 

o 

lA 

O 

lA 

NO 

r>. 

s*  © 

O 

rH 

CA 

CM 

• 

• 

• 

• 

• 

• 

« 

• 

1 CA 

o 

M 

CM 

CM 

rH 

NO 

m 

rH 

m 

CM 

O 

rH  r- 

>4 

CA 

CA 

CA 

rH 

rH 

O CA 
lA 

CA 

> 

A. 


O 

4» 

X 

r4 

rH 

X 

r-4  — 

4) 

« 

« 

«/j  co  » a> 

« 

a n*-4 

A. 

X 

'TS 

A. 

5 5 t*  Dc 

* 

* 

0 

X 

X 

V 

3 

*0 

H 

. 

ao 

25 

o 

z 

« 

O 

O 


£ 


o 

© 

o 

CA 

CM 

*n 

o 

o 

O 

•A 

rH 

© 

© 

*4 

m 

in 

O 

• 

• 

* 

• 

• 

• 

• 

• 

» 

CM 

© 

NO 

as 

O' 

00 

o 

NO 

o 

NO 

© 

CM 

o 

n4 

rH 

O' 

NO 

•= 

CA 

CM 

CM 

Mf 

>4 

CA 

CM 

CM 

CA 

*4 

o 

M 

CA 

rH 

rH 

rH 

rH 

m 

rH 

rH 

i 

<r 

1 

o 

o 

o 

♦ 

o 

CM 

CM 

O' 

ao 

rH 

o 

UA 

© 

iA 

r>. 

1 

o 

o 

CA 

rH 

rH 

NO 

• 

■ 

• 

• 

• 

*■ 

• 

Pm 

o 

o 

CM 

f". 

f- 

rH 

O 

in 

lA 

<4 

tA 

CM 

rH 

in 

pk 

CM 

© 

o 

O' 

O 

r-i 

CA 

CM 

rH 

rH 

CM 

CA 

p 

rH 

o 

_* 

rH 

rH 

o 

rH 

CM 

j 

« 

* 


TABLE  3-75  KNOT  HYDROFOIL  CONCEPTS  - DESIGN  PARAMETERS  AND  CHARACTERISTICS,  COlTT 
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assumed  to  b*  repvesentat' vj  of  waterjet  systems  having  hydrofoil 
?^rut/pod  ram-type  inlata.  . lor  the  aecondary  aystam,  assumed  to  be 
of  the  flush  type,  la  taken  to  be  .5.  Both  primary  and  secondary 
systems  are  asau>so4  to  have  a v axlmum  pump  efficiency  of  90X.  This 
leads  to  an  overall  propul'  We  coefficient  of  .482  at  maximum  speed, 
and  an  overcil  propulsive  coefficient  of  .527  at  the  hullbome  endurance 
speed . 


The  estimated  cha  .acterlscics  Include,  in  addition  to  those 
pres-'  t«-  in  'if  e 1,  th«  shaft  horsepover  required  by  the  pump  at 
take-off  speed „ V V ). 

Kt 


Design  paru-  iters  and  resultant  characteristics  for  the  50,  75, 

and  100  k.v.  SES  ship  concepts  are  shown  in  Tables  4,  5,  and  6, 

respectively.  For  » te  "transport"  mission,  i.e.,  VK  ■ VK  * VK  * 

d max  *|e 

a range  of  500  miles  at  is  selected.  For  the  patrol  mission 

€ 

the  endurance  requirements  are  specified  as  2500  miles  at  12  knots  for 

the  50  knot  SES  concepts  and  1500  miles  at  18  knots  for  the  75  and  100 

* 

knot  SES  concepts  and  a secondary  propulsion  system  is  assumed.  For 
all  cases,  the  number  of  accommodations  Nc  and  the  dry  stores  duration 
D are  assumed  to  be  12  and  7,  respectively. 

O 


The  gross  weights  considered  for  all  SES  concepts  are  again  taken 

as  50,  100,  150,  and  200  tons.  Cushion  length-beam,  2./b,  ratios  for  the 

SES  transport  concepts  are  selected  to  be  2.  However,  for  the  SES  patrol 

concepts,  the  Z/b  is  selected  to  t.*  4 except  for  the  100  knot,  200  too 

concept.  The  loss  coefficient  K^j  is  taken  to  be  .5  for  both  propulsion 

systems  and  the  maximum  pump  efficiency  np  for  both  systems  is  assumed 

max 

to  be  .9.  This  gives  an  overall  propulsive  coefficient  at  maximum  speed 
*527,  which  may  be  somewhat  optimistic.  At  the  pre-hump  endurance 
speeds,  the  overall  propulsive  coefficient  is  also  .527. 

* 

A small  diesel  installation  is  assumed,  but  not  identified. 
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TABLE  5-75  KNOT  SES  CONCEPTS  - DESICN  PARAMETERS  AND  CHARACTERISTICS,  COH'T 


PM 

nO 


PM 

v© 

PM 


WN 

& 


sf 

Pn 


O 

PM 


o 

o 


m 

fN 


o 

NO 


o 

PM 


in 

PM 


NO 

4 


00 

PI 

00 


m 

o 

PM 

PM 


o 

4 


NO 

m 


m 

rv 


m 

PM 


U*N 

4 


*o 

PM 


o 

<n 


o 

O' 


pv 

<r 

O 

in 

o 

o 

m 

nO 

on 

o 

in 

r** 

on 

o 

PM 

O 

4 

4 

ST 


«“M  O 
NO 


0O  in 
• • 
ON  rH 


O 

in 


m m 

pm  m 


4 

m 

- SJ 


O M3 

Mk>  NO 

>4  >4 

"7  >4 


NO 

OD 

o 

m 

a 

o 

in 

P- 

in 

rn 

NO 

03 

o 

PM 

NO 

o 

PM 

«-4 

00 

NO 

1-1 

NO 

PM 

PM 

PM 

-4 

On 

NO 


O 

<n 


O 

o> 


m 

'4 


O 

m 


m 

pm 


m 

oo 


CN 

in 

NO 

r*- 

NO 

r- 

CM 

ON 

H 

O 

o 

o 

O 

vO 

o 

03 

m 

m 

m 

o 

*4 

CO 

03 

o 

ri 

r-M 

•H 

in 

PM 


ON 

*4 


CO 

oo 

On 

4 

o 

o 

*-H 

>4 

4 

PM 

*n 

in 

sO 

PM 

<r 

oo 


o 

m 


m 

PM 


in 

oo 


PM 

o 


PM 

NO 


¥ 8 


O 
O >4 
m pm 


PM 

<0 

r o 
" o 


O NO 
m w-4 


o 

ON 

On  « 


I 

Cfc. 


P: 


O 

o 

Pm 


O 


O 

On  ■* 

ON  c 
I o 

U-  00 
Ou  NO 


PM 

NO  o 
X o 
• o 

H NO 
O -H 


in 


o 

ON 

<*  * 

1 o 

£ § 
Cu 


e 3 


in 

pp  o 
o 

cl 

I-  vO 


PM 


X 

V) 


0 <U 
3 3 


CO 

O 

c» 


« 

u- 

0) 

o 

c 

o 

•H 

o 

H 

X» 

«-M 

Qi 

U 

“O 

3 

t-M 

o 

ii 

£ 

H 

X 

a* 

X 

Cl 

< 

« 

to 

w3 

to 

o 

"1 

1 


32 


O MM 

z o 


10000  14000  15000  5500  10000  14000 


TABLE  6 - ICO  KNOT  SES  CONCEPTS  - DESIGN  PARAMETERS  AND  CHARACTERISTICS 


TABLE  6 - 100  KNOT  SES  CONCEPTS  - DESIGN  PARAKETERS  AND  CHARACTERISTICS .CON >T 
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The  estimated  characteristics  presented  In  Tables  A,  5,  and  6 are 
the  same  as  those  given  in  Table  1,  except  that  the  shaft  horsepower 
required  by  the  pump  at  primary  hump  speed  P(V^)  Is  also  Included. 

RECOMMENDED  VEHICLE  CONCEPTS 


For  each  of  the  3 maximum  smooth  water  speeds  considered  - 50, 

75,  and  ICO  knots  - and  for  each  of  the  2 mission  type  profiles 
considered  - transport  and  patrol  - a vehicle  concept  is  recommended 
for  use  as  reference  craft  for  waterjet  propulsion  systems.  They  are 
described  below. 

50  Knot  Vehicle  Concepts 

The  payload  weight  characteristics  for  all  the  50  knot  vehicle 
concepts,  given  in  Tables  1,  2,  and  4 are  considered  acceptable.  In 
fact,  the  payload  weight  fractions  for  all  these  concepts 

are  between  .08  and  .34.  Since  the  transport  mission  requires  a continuous 
sustained  speed  of  50  knots,  the  hydrofoil  concept  is  recommended  for 
this  mission  profile  as  its  performance  in  a seaway  is  superior  to 
the  other  two  type  concepts.  However,  for  the  patrol  mission  profile, 
the  planing  craft  concept  is  recommended  ns  it  vill  be  cruising  at  a 
moderate  speed  (25  knots  or  less)  most  of  the  time.  Moreover,  it 
represents  the  least  complex  concept  in  terms  of  vehicle  technology 
and  cost.  Table  7 presents  the  recommended  vehicle  concepts. 

75  Knot  Vehicle  Concepts 

The  75  knot  hydrofoil  and  SES  vehicle  concepts  are  presented  in 
Tables  and  5.  Insufficient  payload  capability  is  indicated  for  the 
50  ton  transport  vehicle  concepts  and  the  50  ton  SES  patrol  concept,  i.e., 
the  payload  weight  ie  less  than  5%  of  the  gross  weight.  For  both  missions, 
the  hydrofoil  of  larger  gross  weight  is  again  recomm* nded,  even  though 
its  payload  capability  appears  to  be  somewhat  less  than  that  of  the  SES. 
The  selections  are  presented  in  Table  7. 
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TABLE  7:  RECOMMENDED  VEHICLE  CONCEPTS  FOR  WATERJET 

PROPULSION  DEVELOPMENT  PROGRAM 


VR  50  50 

max 

Mission  Profile  Transport  Patrol 

Recommended  Concept  Hydrofoil  Planing 

Craft 

Grose  Weight,  tons  50-200  50-200 


75  75  100 

Transport  Patrol  Transport 

Hydrofoil  Hydrofoil  SES 

100-200  100-200  > 200 


100 

Patrol 

SES 

> 200 
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100  Knot  Vehicle  Concepts 

In  this  case,  the  only  concept  considered  was  the  SES.  The 
characteristics  obtained  for  the  50  and  100  ton  SES  concepts  (Table  6) 
Indicate  that  these  craft  would  have  no  fuel  or  payload  capability. 
Furthermore,  the  payload  weight  for  the  150  and  200  ton  SES  concepts 
Is  less  than  5%  of  Its  gross  weight.  Thus,  all  of  these  concepts 
are  considered  infeasible.  Viable  SES  concepts  for  these  missions 
would  have  to  be  larger  in  gross  weight.  Table  7 presents  the 
recommendat Ions . 
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APPENDIX  A:  Computerized  Method  for  Waterjet 

Propelled  SES  Vehicle  Concepts 


The  computer  program  "RAVSES"  utilizes  the  method  derived  herein 
to  estimate  the  overall  propulsive  coefficient  and  the  drag  of  an  SES 
vehicle  concept  over  its  operating  speed  range,  and  preliminary  vehicle 
and  component  weights.  The  program  operates  in  either  of  two  modes; 
one  being  specification  of  payload  weight;  the  second  being  specification 
of  gross  vehicle  weight.  When  payload  is  specified,  gross  weight  Is 
derived  by  iteration  until  the  specified  payload  is  satisfied.  When 
gross  vehicle  weight  is  specified,  the  payload  weight  Is  calculated 
directly.  If  the  resulting  payload  Is  less  than  5 percent  of  the  gross 
weight,  the  gross  weight  is  increased  iteratively  until  the  5 percent 
minimum  payload  is  satisfied.  The  program  will  also  calculate  the 
predicted  drag  of  the  vehicle  over  a specified  speed  range. 

A description  of  the  input  to  program  "RAVSES",  a sample  output, 
and  a program  listing  follow.  The  headings  on  the  output  follow 
the  notation  used  in  the  derivation  of  the  method. 
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INPUT  FOR  "RAVSES 


Card 

Columns 

Format 

FORTRAN 

Designation 

Explanation 

1 

1-80 

A 

TITLE 

Any  identification  to  be 
printed  at  top  of  each  page 
of  output. 

2 

1-10 

FI  0.0 

LOB 

Cushion  length-to-beam  r«.tia 

11-20 

F10.0 

LOBSW 

Sidewall  length-to-beam  ratio. 

21-30 

FI  0.0 

WPL 

Payload  weight,  tons.  If  WPL=0.0 
and  DELO  Is  non-zero,  payload 
weight  Is  calculated  from  DELO 
and  other  weight  requirements. 

31-40 

FI  0.0 

DELO 

Vehicle  gross  weight,  tons.  If 
DELG=0.0  and  WPL  is  non-zero, 
gross  vehicle  weight  is 
calculated  from  WPL  and  other 
weight  requirements  by  an  iterative 
procedure. 

41-50 

F10.0 

VKMAX 

Maximum  vehicle  speed  In  smooth 
water,  knots. 

51-60 

F10.0 

VKD 

Vehicle  design  speed,  knots. 

61-70 

F10.0 

VKE 

Vehicle  endurance  speed,  knots. 

71-80 

F10.0 

R 

Vehicle  range,  miles. 

3 

1-10 

F10.0 

NC 

Number  of  accommodations. 

11-20 

F10.0 

DS 

Dry  stores  duration,  days. 

21-30 

F10.0 

KD 

Waterjet  system  loss  coefficient. 

31-40 

F10.0 

ETAPM 

Maximum  pump  efficiency. 

41-50 

F10.0 

PLOPP 

Lift  horsepower-to-propulslve 
horsepower  ratio. 

51-60 

F10.0 

EPS 

Governs  accuracy  of  gross  weight 
calculation  when  DEI.0=0.0  is 
entered.  If  the  difference 
between  the  initial  gross  weight 
estimate  and  the  calculated  gross 
weight  is  greater  than  EPS  tons, 
another  Iteration  is  made. 
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FORTRAN 

Card 

Columns 

Format 

Designation 

Explanation 

3 

(con't) 

61-70 

F10.0 

SYST 

Use  1 for  a single  gas  turbine. 
Use  2 for  a gas  turbine  primary 
system  plus  a secondary  system 
for  low  speed  operation.  If  th< 
endurance  horsepower  Is  greater 
than  3000,  the  secondary  system 
Is  assumed  to  be  gas  turbine; 
otherwise  It  is  assumed  to  be 
diesel. 

71-80 

F10.0 

PRNT 

Use  0 If  each  iteration  is  not 
to  be  printed.  Use  1 if  each 
Iteration  is  to  be  printed. 

4 

1-5 

15 

NVK 

Number  of  speeds  throuqh  speed 
range  for  which  drag  prediction 
are  to  be  made. 

If  NVK  = 0,  omit  cards  4 and  5 
through  T. 

5 

1-10 

F10.0 

VK(1) 

Speeds  In  knots  for  which  drag 

• 

11-20 

F10.0 

VK(2) 

predictions  are  to  be  made. 

• 

• 

Enter  up  to  8 speeds  per 

T 

• 

* 

VK(NVK) 

card  to  a maximum  of  60 
speeds. 

T+l 

1-5 

15 

MORE 

Use  0 if  another  case  follows. 
Use  -1  If  no  more  imput  follows 
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PROGRAM  RAVSCS ( INPUT .OUTPUT , TAPE5* I NPUT . TAPC6*0UTPUT ) 

DIMENSION  VK (60) 

LOGICAL  NET 

REAL  K0»L»L0B*L08SM*NC»NU 

COMMON  /ORAG/OAOM .OCOW «00W .OSAOM.OSKOM .OSVO«*OSWO«»OMOW» N 
COMMON  /PHYS/GPAV.NUtPI.RHO.RHOA 

COMMON  /PWR/ALFA.£NG«ETA.£TAPM*KO.NOEN.PLOPP»PIO*TPMC»PM1*PPOW» 
*PVKt oPVKHopVKH 

COMMON  /S£S/BSW,OS.FRMUMP.FRMAX,L,LOB«LOBSW,NC.PQtPOL*R.SYST, 
♦TITLE  <8> .VKO.VKE.VKMAX 

COMMON  /MGHT/OELO.DELl »WAUX,WCES»P£ »WF . MN.VPLvWS. vsl 

MAXIT-S 

OOPL*2.0 

GRAV*32.174 

NU*1 .281 7E_5 

PI ■3.1*15926535898 

RMO- 1.9905 

RHOA*0. 00238 

01*2.0 

02-1 .6878*1 ,6878/GRAV 
03*02/2240 . •• ( l ./3. ) 

1 READ  (5,500)  (TITLF(I) .1=1 .8) 

SOO  FORMAT (8 A 1 0 ) 

RE  AO  (5,502)  LOB.LOPSW»*<PL,DELO.VKMAX,VKD.VKE*R.NC*DS.KD»ETAPM. 
♦PLOPP.EPS.SYST.PRNT 
502  FORMAT (8F10.01 
WRITE (6,600) 

600  FORMAT (1M1) 

!F(OELO.LE.O.)  MET*. TRUE. 

IF (WPL.LE.O. ) MET*bFALSE. 

IF (MET)  OELO*WPL*DO»L 
FRMUMPbSORT ( 1 ./PI *81/180. •L08«10M) 

POL*01/S0RT(L08) 

N5TEP»0 

2 NSTEP-NSTEPM 
W*DELO*2240. 

FRM?«Q3*Vt<MAX*VKMAX*  (P0L/0EL0/L08)  **  ( 1 ./3.  ) 

FRMAX*SORT (FRM2) 

L*02«VKMAX*VKMAX/FRM2 

PO=POL*L 

BSW-l/LOBSW 

04*1 .6878/SQRT (G«AV*L) 

ALFA*1.0 

IF (SYST .E0.2. ) ALFA=vKO/VKMAX 
CALL  RSACV(VKE*04.PVKE: 

PMC-PVKE 

IF  (FRHUMP.GT . VK0*04)  ALFA  = 1.0 
IF (FRMUMP.lt. FRMAX)  CALL  RSACV (FHHUMP.PVKH) 

alfa.i.o 

CALL  RSACV (FRMAX.PVKM) 

IF (FRHUMP.GE.FRMAX!  PV«H=PVKM 
PM I »PVKM*1 .15 
IF (PVKH.GT.PMI ) PMIsPVKH 
CALL  ENGINE (ENG. PMI.dmC.NOEN) 

IF(MFT)  GO  TO  3 


? 

i 


J 
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wpl-oelo-oeli 

If (NSTfP.EQ.l .OR.PRNT.NE.O.)  call  out<i»^step> 

If  IWPL.QE  .0.05*r»£L01  GO  TO  4 
OELO-1 < 15*0EL I -WPL 
if instep. le.waxit)  wpl-o.o 

If (NSTfP.oT.HAXjTl  *.? 

3 If (PRNT.NC»0. ) CALL  OUT(l.NSTEP) 

If (A8S(0ELl-0EL0>.Le.EfS)  GO  TO  4 
OELO-OELl 

IP(NSTEP.LE.NAXlT)  GO  TO  2 

4 CALL  0UT(?v«5TEP> 

WRITE (6*602) 

602  PORMAT (1H0**V  (KTS)  ORAG (LBS)  fN  ETA  PO»CR 

♦ OWOW  OSAOw  USVOW  OSMOW  OEO« 

♦ OAOW«//» 

ALfA-i.O 

CALL  RSACV(fRMAX.P) 

CALL  LINE  (fRMAX.P.L.GRAV.ETA) 

If  (SYST.EQ.2.  ) ALFA-VKO/VKMAa 
CALL  RSACVIVKD-04.P) 

CALL  LINE  <VKD«04.P,L«GRAV.ETA> 

CALL  RSACV(VKE«04.Pj 

CALL  LINE  (VKE«»04.P.L«GftAV,ETA> 

RE AO  (5.S04)  NV* 

504  fORMAT ( IS) 

If(NVK)  99.1.5 

5 RE  AO  (5.502)  < VK  < I ) ,I«1,NVK) 

DO  6 1=1 .NVK 

PR- VK  ( I ) *04 
ALFA-1. 0 

If (SYST.E0.2. . ANO.VK ( I ) .LE. VKO)  ALP A-VKD/VKMAX 
CALL  RSACV(fR.P) 

CAI.L  LINE  (fR.P.L.GRAV.ETA) 

6 CONTINUE 

READ  (5.504)  MORE 
If(MORF)  49.1.1 
99  STOP 
ENO 


DOW 

OSKOW 
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f SUBROUTINE  R$AC V(FR,P) 

V DIMENSION  TABFR(42) ,TABR<42) 

REAL  KO.L.LOB.LOBSW.NC.NU 

COMMON  /DRAG/OAOW.DFOK  .OOm.DSAOw .USKOM .OSVOM.DSNOV.OmOw. W 
COMMON  /PHVS/GRAV.NU.PI .RmO.RHOA 

COMMON  /PWR/ALFA*ENG.ETA.tTAPM.«O.NOEN,PLOPP.PLO«iPMC.PMl.PPOy. 

♦ PVKf «PVKH«PVKM 

COMMON  /SES/BSW.DS.FMHUMP.FfiMAX.L. LOB. LOBS*. NC*PO»POL*R.SYST, 

♦ TITLE  <B>  «VKO«VKF* VKMAX 

OATA  TA0FR/O , * . 1 * . 1 5 * . 2 « .225  * .25  * .27  * • 2fl * .29* • 3 * .31  * • 3 1 9* 

♦•325. .33* .34 t . 35 *<36 ..373*. 38*. 39 ..4*. 425 *.45. .5* 

♦ .55  ♦ .0 « .65  *.67S*.7«.75*.9«.85*.9*1»*1»1«1.2* 

♦I .3* I .4*1 .5*1 .A* 1 .7* 1 . 0/ 

DATA  TarR/O.  * • 006*  »0  t 7 « .0“*  .677  « • 14.  • I 7 . .2*  .4.  .51  * .603*  »t>2* 

♦ .6 1 ? « .582 .*48fl*.41*.38*.37*.372*.4*.44..6..84.1.22. 

♦1.435* 1.538* 1.588*1. A. 1.597. 1.578 *1.542* 1.49* 1.44 <1.36, 1.3* 1.252* 
♦1.218*1. 188*1. 162. 1.14.1.125*1.111/ 

‘ EUt E»*0. 5772156649 

, EOEP*(SORT( 1 .-«n*2.*(K0*SURT (KO> )/(FR/ALFA/F«MA*>*»2>-l . » 

t ♦ • (fr/alfa/frmax ) ••?/ (ko*sqr  t <kd) ) 

ETa-EOEP*ETaPM* ( l.-SOMT  <KO) / 1 1 ..SORT ( KD  > > * ( 1 * -FR/ALr  A/FRMa* > *«2 ) 

IF (FR.GE*  0 • . AND.FR.LT . 0 *5*FRHUMP ) 

18AVEPAR» 16. •LOfl'*FR*FB/ (PI ®FRMUMP**4 • * ( I .♦ I .A*LOR**.25>  **2> 
2«(3.-EULER*ALOCi<  ( (l.*l.6*L0H**.25>*FRKUMP)4*2/(16.4L0B)  ) > 
IF(FR.GE.0.5*FRHUmP.AnO.Fk.lE.O.B*LOB**.25*FWHUmP> 

1 WAVEPAR=E XP ( -SORT (LOB)*FR«*0.25)*SIN(PI •FRNUMp**2  / (2.*FR««2) ) »«2 
2*LO0/(PI#(FR*.R»LOR*».2S«FRHUMP)»*2) 

3* (3 .-EULER* ALOG ( (FR. ,fl«LOH»*.2S*FRNuMP ) **2/ (4 .*LOB> ) ) 
IF(FR.GT.0.8«LOR**.?5*FRMUMP) 

lWAVEPARaExP (-SORT (LOP) •FR**0.25) *SIN  f PI *FRHIIMP**2/ (2,*rR«*2) ) ••Z 
2*LOB/(4.«PI*FR«*2)  • (3.-LULER  ♦ ALOG (FR**2/LOP > > 
DWO8«4.»WAVEPAR«P0L/RN0/GRAV 
VFS*FR*SO»T (GRAV«L) 

08.RH0»VFS*VFS/2. 

IF(FR.lT.FRHUMP)  HAaPO«(l.-0.8*(FR*FR/FRMUMP/FRhUMP) ) /RHU/GRAV 
; IF  (FR.GE. FRHIIMP)  MA=0.2#Pu/RMn/G«AV 

RE-VFS»L/NU 
CF«FRICT(1 *RE> 

f DSAOW*4.0*08*CF*L»HA/W 

! 0SVO8»OW*CF*LOR«O80W/P0 

oi«b.o*rho#grav/pi 

IF (FR.LE* 1 *8)  CALL  niSCOT(FR*FR*TABFR*TABR*TABR»-iaO*42.0.Rl) 
IF1FR.GT.1 .8)  Rl*6.125«ALUGtFR>/FR/FR 
OSWOWsQl  ®R1  *BSW#0SW«HA*HA/IK/L 

qa*rhoa»vfs*vfs/2. 

0EON»0.2#OA/L0R/P0 
OSKON-O . 1 •CF*OW/PO 
COA»0.393E~4*0.02*CF 
DAOW“COA#OW/P0 

DOW«DWON*OSA08*OSVOW*OSW08*OE08*OSKOta*DA08 
PPON*OON»VFS/ETA/550 . 

PLOW*PLOPP#PPOW 

P«(PPOW*PLOW)*W 

RETURN 

END 
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SUBROUTINE  ENGINE (ENG* PM  I .Rmc.N) 

DIMENSION  T ARENA ( 3* I I > »NE 111) .NMaM  10) .PMCNlll) .PMINU  1 > 

COMMON  /WGHT/OElO»nFLl »WAUX.WCES*WEtWF.WM«WPL.WS»iSL 

DATA  TaRENG/IOH  G NOME. 8?5. *970. *10H  LM100U150. . 1220. . 

♦ 1 OH  TE1AR  * 1 375. • 1550* • 10H  Te25a*2200.»2500.*1Qh  TF3S 

♦2?S0.»3150.» I0H  S01-K1A.1A00..3730..IOH  <jTpr-v*0» 5000* .5600. t 

♦ 1 OH  501-M62.7000.. "000. »10h  9I2-C1 . 1 3000. • 15000. • 

♦ 1 Oh  lm2500*?a050.»?$600..10h  ET9a-2»36000. .A2500./ 

OAT*  NMAX/6*(?) .3«(3> *5/ 

no  i i-i.u 

NC  ( I ) -POUND  HPMI  /T ABENG  < 3. 1 ) ♦ .5)  »1.) 

IF(Nrm.lT.2>  NE ( I ) «? 

IE  <T ABENG (2. I ) *NE ( I ) .LT.PMC)  NE ( 1 > -NE ( I ) » 1 
PMCN ( I ) -t ABENG (?» I ) -NE ( t ) 

1 PMIN(I)-TABENG(3.I)*NE(I> 

IMIN-Il 

CALL  WEIGHT  <NE<U)  .PMCN  ( 1 I ) .PMIN(II)  > 

OELMrOELt 
00  ? 1-1. 10 

IE(NE(I).GT.NMAX(In  GO  TO  2 
CALL  WEIGHT INE ( I) »PMCN( I ) .PMINU) ) 

IE (DELI .GT.OELM)  GO  TO  2 
IMIN*I 
DELM-OEL1 
? CONTINUE 
N-NE ( Imin> 
eng-tabfngu  .IMTN) 

PMCaT  ABENG ( 2. I MIN ) **N 
PM  t * T ABENG  ( 3 • I M I H ) **N 
CALL  WEIGHT (N, PMC. PMI) 

RETURN 

FNO 


SUBROUTINE  WEIGHT <NE»PMCNiPHIN> 

REAL  L'LO0»LO(tSW«NC 

COMMON  /PWR/ALrA«ENfl»ETA«ETAPM*KD»NOENtPLOpp»PLOW*PMC*PHl*PPOW* 
♦PV*E»PVKM*PVKM 

COMMON  /SCS/8SW*0S*rRHUHP«*'RMAX«L iLOBeLOBSw»NC«PUIPOL«R*SYST* 
♦TITLE* A) fVKDtVKE »VKMAX 

COMMON  /WOHT/OELO*OEL  1 » WAUX  , WCES » WE  t WP  ♦ WM  * Wpl  ♦ WS  ♦ R!»L 
lr <RVKE.OT*3OO0**OR,SYST.EQ«l.)  go  TO  | 

SPC«r0.4 

WN*0.001562*PNfN 
80  TO  2 

I NbROUNO ( ‘ °VKC/PHCN»NE ♦ • 5 ) « 1 . ) 

PEN*RVXE/N 

SPC«1.9»PCN«#(-0.1S) 

WN«0.001339«PMIN 

t WS-OELO/SORT (POL ) • <OELO«« ( -1 «/3. ) *0 . 3/L08»» (1 »/3«T> 
WSL«0.05»OELO 
WC»O.03#r>CLO 
WAUXaO.lPDELO 
WCCS-0 . 1 *NC*0 . 0 15*NC*0S 

Wr*OELO»(l.-EXP(-I.«SrC*H»PVKE/(22AO.«OELOavKE) ) » 

0EL1«WPL»WS*WSL»WM«WE«WAUX»WCES*MF 

RETURN 

END 
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SUBROUTINE  OUTUPAGE.NSTEP) 

REAL  KO»L»LOfl»LOBSW.NC 

COMMON  /PWR/ AIF  A «ENG  »E  T A • t T APM.Kl)  »NOEN  .PLOPP*  PLOW  I PMC  »PM  I *PPOW  • 
♦PVKF.PVKH.PVKM 

COMMON  /SES/flSW  »OS»FRhUMP»FRmaX ,L  *LOB«LOBSW»Nc»POJ POL*R*SY$T  , 

♦title (A) .vkd.vke* vkmax 

COMMON  /WGNT/DELO»0fLl***UX.WCES.WE»WF,WM»WPL.WS»»5L 
GO  TO  (1.2)  IPAGE 

1 WRITE (6.600)  (TITLE(I) *l»l.B) 

600  EORMAT(////28X.ftA10! 

GO  TO  3 

2 WRITFI6.602)  ( T I TLF ( I ) ,1*1.8) 

602  rORMATdHl  .2?X,flAlO) 

3 WRITE(6,604)  NSTfP 

604  FORMAT (1H0.25X,«RAVENSCR0FT  METHOO  FUP  WATERJET  PROPELLED  SES  CONC 
♦EPTUAL  OE5IGN*. 10X, •ITERATION  *ill) 

WRITE (6,606)  VKMAX , VKO. VKE ,R .NC.OS »LOB .POL . PO.L 
606  FORMAT (1H0.«  VKMAX  VKO  VKE  RANGE  ACCOM. 

♦ DRY  ST.  L/R  P/L  PO  LENGTH«/1X,7F10.0.P10.3. 

♦F10.0.F10.' ) 

FRE«VKE»FRMAX/VKMAX 

WRITE (6.608)  NOFN.FNG.PMC.PMI .FRfc ,PV«E .FRMAX .PVKM* FRHUMP 
608  FORMAT (1HO.*  NO.  ENG.  MOOEL  PMC  PMI  ERE 

♦ PVKE  FRMAX  PVKMAX  FRHUMP  PVKM  KD  ETA 

♦ ETAPMAX*/! 1 l,A10.2E10.0.3(F10.2.riO.O) ) 

IF (FRHUMP.LE. FRMAX)  WR I TE (6 .6 1 0 ) PVKH 

610  FORMAT (lH*.R0X»FlO,0) 

IF (FRHUMP.GT. FRMAX)  WRITE (6,612) 

612  FOPMAT(1H..90X,*  — — •)  • 

WRITE<6.614)  KO « E T A « ptapm 
614  FORMAT (1H..100X.3F10.2I 

WRI TE (6 .61 6)  OELO.nELl .WS.WSL.WM.WE.WAUX.WCES.WF.RPL 
616  FORMAT ( 1 HO » * OELO  OEL1  WS  WSL  WM 

« WE  WAUX  WCES  WF  WPL*/1X.2F10,0.8F10.1) 

RETURN 
END 
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SUBROUTINE  LINF<PR.P.L»G*F.TA) 
real  L 

common  /ORACi/nAOW.OEOW»D<>»»OSAOW.OSKOW.nSVOW»DSWOW»OWOW.W 

Vk*FR*SQRT  (G#L)  /I  .6B70 

o*oow*w 

WRITE  (6 #600  Vk#O»PR#ETa*P#DOW#0*(0W»0SaOW#0SVOW*0SW0W»0EOW#DSKOW* 

♦OAOW 

600  format <lX.E6.2.1P£ll.3.0P2F7.2*lP9tll.3) 

RETURN 

END 


SUBROUTINE  OISCOT  (XA.ZA.TABX.TA0Y.TAftZ.NC, NY. NZ.ANS) 
DIMENSION  TA0X(  10)  .TA8T  ( U<#0)  • lASZ  ( 14A)  ,NPX  (8)  #NPT  (0)  »YY(8) 

OISCOTOl 

CALI.  UNS  (NC#  IA,  IOX#  I0Z»  IMS) 

DISCOT03 

IF  (NZ-1)  5#5#10 

DISC0T04 

5 

CALL  DISSER  ( XA . TARX# 1 *NY. IOX.NN) 

OISCOT05 

NNN« IDX ♦ 1 

DISCOT06 

CALL  LAGRAN  (XA.TABX(NN) .TA8YINN) .NNN.aNS) 

OISCOT07 

GOTO  70 

D I SCOTOfl 

10 

zaro>za 

DISC0T09 

IP1X*I0X*1 

DI SCOT  1 0 

IP1Z*IDZ*1 

DISCOT  1 1 

IF  (IA)  15.25.15 

OISCOT 12 

IS 

IF  (ZARG-TABZ(NZ) ) 25.25,20 

OISCOT 1 3 

20 

ZAPG«TA0Z(NZ) 

D I SCOT  14 

?5 

CALL  DISSER  (ZARG.TABZ.l .NZ, IOZ.NPZ) 

DISC0T15 

nx«ny/nz 

OISCOT  16 

NPZL*NPZ ♦ IOZ 

DI SCOT  1 7 

1*1 

DISCOTlfl 

IF  (IMS)  30.10.40 

DI SCOT  19 

30 

CALL  DISSER  (XA.TARX.l.NX.IOX.NPX) 

DISCOT20 

DO  35  JJaNPZ.NPZL 

DISC0T21 

NPY ( I) * ( JJ-1 ) #NX  *NPX ( 1 > 

DISC0TZ2 

NPX ( I ) *NPX ( 1 ) 

DISCOT23 

35 

1*1*1 

DISC0T24 

GOTO  50 

DISC0T25 

40 

DO  45  JJ=NPZ.NPZL 

DISC0T26 

IS*(JJ-1 ) *NX*  1 

0ISC0T27 

CALL  DISSER  (XA.TABX.IS.NX. IDX. NPX ( I ) ) 

OISCOT20 

NPY(I) *NPX ( I ) 

DISCOT29 

45 

1*1*1 

OISCOT30 

50 

00  S5  I*1,IP1Z 

DI SCOT 31 

NLOC«NPX (I) 

DISCOT32 

NLOCY*NPY(I> 

D 1 SCOT  33 

55 

CALL  LAGRAN  (XA.TAPX (NLOC) .TA0Y (NLOCY) . IP1X.YY ( I ! ) 

01  SCOT 34 

CALL  LAGRAN  ( ZARG . TA0Z ( NPZ ) . Y Y ♦ IP  1 Z . ANS ) 

DISCOT35 

70 

RETURN 

DISC0T36 

END 

0ISC0T37 

48 


I 


SUBROUTINE  LAORAn  ( X A * X . Y «n.  AnS  ) 
DIMENSION  XUO)  .¥(1440 

0ISC0TS4 

SUM»0.0 

0ISC0T56 

no  3 I»1,N 

OISCOT57 

p«no*y ( I ) 

01SC0TS8 

no  ? jsi.M 

01SC0T59 

A«X(D-X(.n 

OISCOT60 

If  (A)  1.2.1 

0ISC0T61 

1 

R«<XA-X(J» )/A 

0ISC0T62 

ppnn*PRon*R 

01SC0T63 

? 

CONTINUE 

01  SCOT 64 

3 

sum=sum*proo 

0ISC0T65 

AN5*SUm 

OISCOT66 

return 

0ISC0T67 

Eno 

0ISC0T68 

SUBROUTINE  OISSER  < XA , TAB. 1 .NX. IU.NPX) 
DIMENSION  TAB(IO) 

DISCOT69 

NPTsin.l 

DISCOT71 

NPRaNPT/? 

OISCOT72 

npu»npt-npb 

OISCOT73 

IE  <NX-NPT>  10*5.10 

DISCOT74 

S 

NPX*  I 

OISCOT75 

RETURN 

DISCOT76 

10 

NLOW* 1 »NPfl 

OISCOT77 

NUPP«I*NX-<NPU«! > 

OISCOT78 

DO  15  TI^NLOW.NUPP 

DI SCOT  79 

NLOC-1 I 

DISCOTBO 

IE  (TAP(IT)-XA)  15.20.20 

d i Scots i 

15 

CONTINUE 

discotsz 

NPX*NUPP-NPB»1 

DISCOT83 

return 

DISC0T84 

20 

NL«NLOC-NPB 

0ISC0T85 

NU«NL‘ID 

0ISCOTS6 

DO  25  JJ=N|_»NU 

DI SCOT  87 

N0I5«JJ 

0ISCOT8B 

IE  <TAR«JJ>-TAB< JJ*1) ) 25.30*25 

DISCOT89 

?5 

CONTINUE 

DISCOT90 

NPX«NL 

0ISC0T91 

return 

0ISC0T92 

30 

IE  (TA8(NDIS>-XA)  40.35.35 

DISCOT93 

35 

NPX*NDIS-IC 

DISC0T94 

return 

DI SCOT  95 

40 

NPX*ND I 5. 1 

D1SC0I96 

RETURN 

END 

DISC0T97 

49 


5 

J 


SUBROUTINE  UNS  (IC»IA*IOX.IDZ.lMS>  OISCOT30 


if  do  s«5.io 

OISCOT39 

s 

IMS-1 

OISCOT40 

NCtf-IC 

DISC0T41 

GOTO  15 

0ISC0T42 

10 

IMS-0 

0ISC0T43 

NC-IC 

0ISC0T44 

IS 

I-  (NC-100)  20*25.25 

DISCOT45 

20 

IA-0 

0ISC0T46 

GOTO  30 

01 SCOTS? 

*5 

IA-1 

DISC0T4B 

NC -NC-100 

0ISC0T49 

30 

IOX-NC/10 

OISCOT50 

IOZ-NC-IOX-IO 

DISC0T51 

RETURN 

0ISC0T52 

ENO 

DISCOT53 

FUNCTION  FRICT(IFRICT.RE) 

FRI 

001 

GO  TO  (4.2)  IF9ICT 

FRI 

002 

2 

FRICT-.075/1ALOG10 (Rf/100. 1 ) **2  SGO  TO  99 

FRI 

003 

4 

X-l  ./(3.46-ALOG10  (RE) -5. 6) —2 

FRI 

004 

6 

FPTCT*(.242/ALOG10 (RF-X) > t-2 

FRI 

005 

ir<AB5(FRICT-X» .LT.5.E-07)  GO  TO  99 

FRI 

006 

Xa.  <FRICT*X ) /2» 

FRI 

007 

GO  TO  6 

FRI 

008 

99 

RETURN 

FRI 

009 

ENO 

FRI 

010 

FUNCTION  ROUND ( X • V ) 

RNO 

(*" 

C 

X IS  VARIABLE  TO  be  rounded 

RND 

002 

c 

V IS  The  value  TO  WHICH  x is  to  be  rounded 

RNO 

003 

XliX/V  STX-NX-X  SDX*X-TX  STdX*NX*2.0«DX  SROUND-  <TX*TDA)*V 

RNO 

004 

return 

RND 

005 

END 

RND 

006 

50 


